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usando OpenFoam 


Abstract 

Biomass is an important renewable energy source that has great potential as a substitute 
for fossil fuels in the short and medium-term, which has led to the development of 
various methods for its energy conversion, of which combustion is the most widely 
used. This process has several environmental advantages compared to traditional 
energy sources, however, there is still a long way to go in terms of process efficiency and 
emission reduction. In this context, CFD computational models are a powerful tool that 
enables the study and improvement the performance of combustion systems in a safe, 
fast and economical way, compared to experimental studies. OpenFoam is one of the 
most important CFD software currently available, yet there are few works that use it to 
simulate the combustion of solid biomass. In this work, an application of the software 
in modeling of a biomass boiler fueled by grape marc is reported, this model allows for 
the prediction of important parameters like CO, H2O and velocity fields in a 2D domain. 
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S) ® © Resumen 


La biomasa es una importante fuente de energia renovable que tiene un gran potencial 
como sustituta de los combustibles fdsiles en el corto y mediano plazo, lo cual ha llevado 
al desarrollo de varios métodos para su conversion energética, de los cuales la combustién 
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Received: la eficiencia del proceso y la reduccién de emisiones. En este contexto, los modelos 
WEA GE computacionales CFD son una herramienta poderosa que permite estudiar y mejorar el 
Aceptado / desempeno de sistemas de combustidn de una manera segura, rapida y econdmica, en 
comparacion con estudios experimentales. OpenFoam es uno de los softwares CFD mas 
importantes disponibles en la actualidad, no obstante, hay pocos trabajos que lo utilicen 
Publicado en linea / para simular la combustidn de biomasa sdlida. En este trabajo se reporta una aplicacién 
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orujo de uva, este modelo permite predecir parametros importantes como el campo de 
velocidades y las emisiones de CO y H20 en un dominio 2D. 
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INTRODUCTION 


To overcome the serious environmental problems of global warming and poor air quality, 


it is necessary to develop clean energy sources and reduce society's energy consumpti 


on. 


Studies have been carried out around the world in order to find and improve renewable 
energy sources that can replace fossil fuels. Biomass is one of the most important; it is 
considered the most promising energy source to mitigate greenhouse emissions in the short 


and medium-term [1]. The use of biomass-derived power has increased during the last 
g the fourth largest energy source worldwide, after coal, oil, and natural 
gas [3] [4] [5]. However, there remain challenges associated with optimizing and applyi 
s which make biomass energy conversion seemingly too inefficien 
\y and commercially viable alternative to traditional finite fuels [2]. 


decades [2], becomin 


biomass technologie 
be an environmental 


few 


ng 
to 


Over time, various processes have been developed for the transformation of biomass into 


energy, framed 


within thermochemical, biochemical and chemical conve 


rsion routes. Of 


these, combustion is the most used as it is the most suitable from the economic point of 


view [6] 
such as slaggin 


necessity in terms of resea 
s, as the increase 
consumption; analysis has shown that a 1% in 
plants causes a 3% reduc 


pollutan 


needed in the 
properties, whi 


of the use of other types of biomass is expected, 
and biogenic waste [9]. These will require the s 


1]. Nevertheless, t 


g, instability of heat load, and low thermal efficiency [7] [8]. 


in efficiency leads to a significant decrease in em 


here remain problems in current biomass combustion furnaces, 


hese present a 


ch to improve the performance of these devices and to reduce 


issions and fuel 


crease in the efficiency of biornass power 
ion in their CO2 emissions [1]. Another area of 


improvement 


research is the biomass characterization in terms of its physico-chemical 
ch is imperative to understanding biomass combustion. Traditionally the 
most commonly used solid is wood but, in order to improve energy utilization, an increase 


inc 
udy of their properties and 


processes to achieve maximum energy performance. 


Computationa 


models are a valuable tool 


uding a diversity of agricul 


ural residues 
transformation 


o improve the performance of biomass 


combustion systems since - the design and operating conditions of these systems can be 


improved with 
studies[10] [11 
(CFD) models i 
Dernbecher et 
is an issue tha 
process, which 
biomass and s 


out incurring the economic and environmental costs of th 
. This has motivated the development of Computational 
n recent decades. A full review of this topic is presen 


requires various submodels for the thermo-chemica 
ufficient computational resources for the simulation 


ed i 


e experimental 
fluid dynamics 
in the works of 


al. [6], Karim and Naser [12], Bhuiyan et al. [1], and Khodaei et al. [2]. This 
is still under development due to the complexity of the combustion 
conversion of the 
9]. OpenFoam is a 


powerful open-source engineering tool, which has been successfully used in numerous 
applications. It has the advantages of being a free and open source, being constantly 


under development, and including solvers for many applications 


13). 
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rs, there are relatively few works that use it to simulate the combustion of solid 


facto 
| biomass. In order to guide researchers interested in the use of OpenFoam for this type of 
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cations, this document provides a 2D simulation of combustion in a boiler fed with 


grape marc employing OpenFoam [14] [15] [16]. This is a simple model that predicts CO 


and H,O profiles of velocity and temperature over the freeboard. 


METHODOLOGY 


high 


he s 


see 


is work, the combustion in a biomass boiler with dimensions 8 m wide and 15 m 


ina 12.5 MWel combined heat and power plant was simulated. The system was 


already analyzed in the work of Cordiner et al. [17] [18], but the novelty of our study lies in 


implicity of the model, which requires relatively few computational resources. The 


furnace operating conditions are presented in table 1. A squared mesh with a grid size of 
approximately 20 was used, the computational grid was done with the blockMesh utility 


figure 1). The characteristics of the biomass used are shown in table 2, which were 


also 


aken from the work of Cordiner et al. [18]. The volatilization kinetic parameters were 


aken from the work of Fiori [19]. 


Table 1. Furnace operating conditions [18]. 


Condition Value 


Biomass flow rate 18,300 kg/h 
Primary air flow rate 53,750 kg/h 
Secondary air flow rate 53,750 kg/h 
Secondary air nozzles diameter 76mm 
Secondary air nozzles number 38(19+19) 
Primary air temperature 170°C 
Secondary air temperature 25°C 


Table 2. Fuel characteristics [18]. 


Property Value 


Volatile Fraction 30% 
@) 80% 
CH, 10% 
Cink. 10% 
Solid Fraction 20% 
Ashes 30% 
C content 52.1% 
H content 5.8% 
N content 0.84% 
S content 0.19% 
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Figure 1. Computational Domain 


The simulation was under the following main assumptions: 


+ The fuel (biomass) is supplied to the combustion chamber following a predefined 
path. 


+ Solid particles are not moved by gas action (drag and gravity forces are neglected) 
but gas phase movement is affected by bed particles. 


» The Equations of continuity, energy and momentum of the gas are solved taking 
into account chemical reactions and interactions with the solid phase. 
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An empirical 2D model implemented in OpenFoam was used, utilizing the 
coalChemistyFoam solver. To model the turbulence/chemistry interaction in the 
combustion process the k-e model and the Partially Stirred Reactor (PaSR) combustion 
model were employed, in the latter, the computational cells are divided into two zones, 
one in which all reactions occur and the other in which no reactions occur so that the 
composition changes due to the mass exchange with the reaction zone. The governing 
equations used in the gas phase are presented below: 


* — Continuity equation: 
200) + 7 - (pgU,) = Sm (1) 


. Energy transpor equationt 


eel + 204K) +V-(pgUgh) + V+ (pgUgk) —V- (Oer;V(h)) =—-Vp + pgUg-g+Sp (2) 


° Momentum equation: 


0 qUy 
Oe sy. (PgUgUg) — V° (tg) — V° (PgRg) = —Vp + pte + Sy 8) 


. Species transport equation: 


Oe V = (pgU,¥;) — V+ (DepV(p+¥3)) = 5S; (4) 


h is the enthalpy, S is the source terms, K is kinetic energy, Rg is Ronald stress term, 
Yi is the mass fraction of the species, and aeff and Deff are the effective thermal and 
mass. The freeboard reaction mechanism used was proposed by Gdémez et al. [20], 
which consists of the reactions summarized in table 3. The system was simulated in a 
computational cluster of the “Universidad Industrial de Santander”, named GUANE, that 
consists in 16 ProLiant SL390s-G7 computing nodes such as shown in Table 4 [21] [22]. 


Table 3. Combustion reactions [21]. 


Homogeneous reactions 


9 
RI CH, + 0, > 6CO + 3H,0 
3 
R2 CH, +, 0, > CO +2H,0 
1 
R3 H, +0, > H,0 
1 
= cO+, 0, >C0, 
RS H,0+CO>CO,+H, 
R6 CO, +H,>H,0+CO 
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Table 4. Se 


Setting Guane 


Node type 
Number of nodes 


Processor Intel 
Processor Model 
Processor by node (#) 
Clock frequency (GHz) 
Core/Processor (#) 
Thread/Core (#) 


GPUS Nvidia 
GPUS Model 
GPUS by nade (#) 


Memory DDR (GB) 
SAS disk (GB) 


Gigabit Ethernet (Gbps) InfiniBand (IB) 


RESULTS 


The furnace was simulated under s 
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tting GUANE Cluster Nodes [21] 

A B Cc 
SL390s SL390s SL390s 
8 3 5 
Xeon Xeon Xeon 
E5645 E5645 E5640 
a 2 2 
240 240 240 
6 6 4 
a 2 2 
Tesla Tesla Tesla 
M2075 M2075 M2075 
8 8 8 
104 104 104 
200 200, 200 
101 101 101 


teady state conditions, CO, H2O and velocity fields 


obtained are reported in figure 2. The results are similar to those obtained by[18], who 
simulated the same type of boiler. Although this model does not simulate NOx emissions, 


its value is that it achieves good resul 


sin asimple way with low computational resources. 


The calculated temperature closer to the outlet section equal to 1300 K, which can be 


compared with the experimental on 
and velocity values at the outlet of t 


e, that is equal to 1,270+15 K [18]. Also, the CO, H2O 
he chamber, evaluated in this study and taken from 


the literature [18], were compared; which can be seen in table 5. The difference between 


results is less than 6%, which shows 
the combustion of solid biomass. Th 


he capacity of the model to predict the behavior of 
erefore, the model developed is useful to study the 


feasibility of the use of new biomass fuels. 
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Figure 2. Simulation results. A) Velocity magnitude field [m/s]. B) CO mass fraction. C) H20 mass 
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Table 5. Comparison between data obtained from the literature and data obtained in this work. 


Parameter Literature [18] 
Velocity [nm/s] 7.50 TOF 
CO mass fraction 0.00 0.00 
H20 mass fraction 0.13 0.10 
CONCLUSION 


A 2D model for biomass furnace is reported in this paper, the model represents the 
evolution during thermal decomposition in the combustion chamber with a relatively low 
computational cost. The model is capable of predicting the CO and H20 concentration 
and velocity, and temperature, profiles evolution over the freeboard. The results present an 
overall satisfactory behavior showing good agreement with literature data. 
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